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In this paper we present the generalized phase separation diagram for quantum plasmas revealing 
the fundamental mechanism called the spontaneous core collapse (SCC) for a white dwarf star core- 
density based on the extended Shukla-Eliasson force for relativistically degenerate matter taking into 
account the most important electron interaction features such as Fermi-Dirac statistical pressure, 
Coulomb, exchange-correlation and quantum tunneling effects. A critical core mass-density value is 
reported for white dwarf stars beyond which the SCC mechanism will be inevitable. It is revealed 
that the SCC critical core-density value is slightly affected by the ion charge state, Z (the core 
composition). Such fundamental effect which is mainly caused by the electron quantum diffraction 
effect can progressively lead the white dwarf star with a mass beyond the Chandrasekhar limit 
through SCC, core fusion and the ultimate supernova stages. 
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INTRODUCTION 



During the past decade, a new dimension of research has initiated due to Bohmian collective interpretation of 
quantum systems [U-EI]. Foundations of quantum hydrodynamic (QHD) formulation, earlier suggested by Madelung 
[J] termed as hydrodynamic formalism for quantum mechanics, was extended by Takabayasi |5j-|9| and others. It has 
been shown that in quantum system of Fermionic particles a new force due to the interaction of particle wavefunction 
may result in important collective plasma dynamic features. Such ne collective feature is called the Bohm quantum 
force 13, 11 1 and is the main ingredients of the quantum hydrodynamics (QHD) formulations. In dense plasmas where 
the inter-electron distances are comparable to or lower than the de Broglie thermal- wavelength A^> = hj (2T7m e kBT) 1 ^ 2 , 
quantum effects such as the Pauli exclusion [l2[ and the electron quantum Bohm force will be inevitable. Quantum 
hydrodynamic theory has been successfully employed to interprets the collective dynamics of electron-hole plasmas 
in a semiconductor [13l4l5j . Variety of applications may be foreseen for quantum plasmas other than semiconductors 
[Tjj ]. For instance, the physical characterization of nanofabricated devices fundamentally rely on the quantum theory 
of hydrodynamics. 

A large number of recent investigations based on QHD model confirm signicant role of quantum diffraction effect 
on linear and nonlinear wave dynamics of degenerate plasmas (l7l - [22l |. It has been shown that the inclusion of 



Bohm force can lead to many critical features [23, l24[. More recently Shukla and Eliasson by investigating the 



screening effects in a quantum plasma have discovered a new attractive force between ions [251 ] which is due to 
interplay between electron degeneracy, exchange-correlation and quantum diffraction effects. The existence of such 
attractive force among ions, if established experimentally, can mark a distinguished milestone in the history of plasma 
physics research with the Bohmian collective interpretation of the quantum particles. Moreover, the Wigner-Poisson 
formulation of hydrodynamics has been recently extended to include the electron spin magnetization effects 26|, [27 1 
or the relativistic electron momentum effects [28(. Using the quantum magnetohydrodynamic (QMHD) theory, it has 
been confirmed that the electron spin can play signicant role in linear and nonlinear wave dynamics of magnetized 
quantum plasmas [IIS]. 

One of the important applications of QHD theory regards the astrophysical applications of degenerate plasmas 



which are ubiquitous in the nature [35H40J • Chandrasekhar [41J, |42| has shown that the compact stars called white 



dwarfs can be effectively modeled as the ideal zero-temperature quantum plasmas. It has also been shown that 
the large external gravitational force on stars can lead to a phenomenon known as the relativistic degenera cy |43l| 
leading to distinct features in equation of state (EoS) of degenerate matter |44j and nonlinear wave dynamics |45| in 
different density regimes. However, the electron diffraction mechanism is supposed to be affected by the relativistic 
plasma effects, since, the relativistic effects are directly coupled to the plasma number density in the relativistically 
degenerate matter. Mettes and Sorg [3] and Carlos and Mahecha 47-50] have extended the fluid dynamics to include 
the generalized relativistic Bohm force by solving the time dependent Dirac or Klein-Gordon equations. More recently, 
a relativistic fluid model based on Wigner function and Dirac equation has been proposed [5l| which includes the 
relativistic corrections to quantum diffraction effect. 

In the following, we consider a zero-temperature relativistically degenerate collisionless electron/ion fluid with non- 
degenerate ions of charge state, Z. The relativistic quantum hydrodynamic fluid model incorporating the relativistic 
electron motion, quantum statistical pressure, interaction effects such as Coulomb, exchange-correlation, etc. and the 



relativistic quantum Bohm force [51|, may be written as 



^ + V-( 7 nu) = 0, 
7(§r + (u-V)u) 

= ^-^vp G + ^v(^f), 

A(j)(r) — Ane(jn — n ), 



where, Pg is the generalized pressure due to relativistic electron degeneracy in addition to other electron interaction 
pressures [52|. Also, 7 = 1 / ^/l — f3 2 is the relativistic factor with /3 = u/c. It should be noted that the additional 
correction to the relativistic Bohm force given in Ref. [5l| is ignored for being relatively small proportional to (3/c. 



Using the concept of generalized effective potential, same as in Ref. 
to harmonic perturbations via the Fourier analysis of Eqs. ([T]) as follows 



531 ] , one may obtain the plasma dielectric response 
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-ami + k • ui = 0, 
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ui = -^jfc 2 kni - Lc 2 kni, ^ 
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It is evident that the relativistic parameter, 7, does not enter the linear response model. In the linearized model, 
Eq. ([2]), the parameter Rq = (jiq/tIc) 1 ^ with n c = 5.9 x 10 29 cm -3 is the Chandrasekhar's relativity parameter, 
a = e 2 /hc ~ 1/137 is the fine-structure constant and (3 = (2ao/5)(2Z/3y / 7r) 2 / 3 . Therefore the static dielectric 
constant is then given as 



D(0,k) = l+iol 



h 2 k± 

Am 2 



Lc 2 



(3) 



where u pe — Aim^e 2 /m e is the electron plasmon frequency. The dielectric around a point-like stationary charge, 
Q = Ze, located at the origin is given by 



Q f exp(ik • r) s 



^ v ' 2tt 2 J fc 2 L>(0,k) 
Decomposition of the integral with known solution given in [25| reads as 



d A k, (4) 



(r) = g[(l + b)e- k - r + (l-b)e- k + r ]. (5) 



with the potential parameters defined as below 

t 2 2 



(6) 



One of the fundamental characteristics of the Shukla-Eliasson potential 25j given by Eq. ([5]) is that for a > 1/4 the 
potential is attractive and repulsive otherwise. The generalized degenerate matter phase separation then is obtained 
by evaluation of the parameter, a. In Fig. 1 we have shown the variation of this parameter for wide density parameter, 
r s = tb/tq range, where, tb = h 2 /m e e 2 ~ 5.3 x 10~ 8 cm is the Bohr-radius and ro = (3no/47r) 1 / 3 is the Wigner-Seitz 
cell radius defined as a sphere containing only a single electron. The relativity parameter is given as Rq ~ 0.014r s . 

As it is observed form Fig. 1, increasing r s (increasing the plasma density) from zero-density first leads the pure 
Coulomb screening (a — +00) to turn over to the Lennard-Jonse-like Shukla-Eliasson attractive potential of form 
4>(r) = (Ze/r) exp[— iJi(fc_)r]{cos[3 , (fc_)r] +6sin[Sy(fc_)r]}, where and 9 refer to the real and imaginary parts of the 
parameter For two different values of r s i(Z) and r S 2(Z) we have a — 1/4. For r s i(Z) < r s < r s i(Z) the plasma 
undergoes the ordinary Debye-Hiickel screening and for r s > r S 2{Z) the Shukla-Eliasson attractive force reappears. 
It is remarked that, the value of r S 2(Z) unlike that of r s \(Z) is not strongly sensitive to the ion charge state, Z.This 
value coincides exactly with the core-density of a typical white dwarf star (10 s < p c (gr/cm 3 ) < 10 13 ). It is remarked 
that r s i = 6.8326 value for iron plasma composition extends to density of p ~ 821.945<?r/cTO 3 which is well beyond 
the zero-pressure metallic iron. This is the indication that the Shukla-Eliasson attractive potential leading to the 
plasma crystallization can be present at small-planet core density ranges. On the other hand, the r S 2 ~ 2354.8 value 
for the same plasma composition corresponds to the mass density of p ~ 3.5 x 10 10 gr / cm 3 which resides in the range 
of a typical white dwarf core. The reappearance of Shukla-Eliasson attractive force in white-dwarf core density has 
some fundamental consequences for stellar evolution, as will be discussed next. 

Figure 2 reveals some fundamental differences of Shukla-Eliassom attraction between nonrelativistic and relativistic 
degeneracy cases. Figures 2(a) and 2(b) compare the potential profiles (normalized binding energy, rB<j>{r)4>/Ze) in 
nonrelativisic and relativistic degeneracy cases against the variation in the ion number density, Z, where f = r/rg. It 
is revealed that, the increase in the number-density leads to increase in the potential depth in both degeneracy regimes 
while this variation is more pronounced in the nonrelativistic degeneracy case. However, as it is remarked comparing 
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FIG. 1. The variation of a against the plasma density parameter, r s — as/ro. The critical value a — 1/4 is indicated with the 
horizontal line. The thickness of curves appropriately reflect the variation in the ion charge value, Z 



Figs. 2(c) and 2(d) the increase in the plasma density causes reverse effects in the two density regimes. In other 
words, the attractive potential strength (binding energy) increases/decreases for the case of relativistic/nonrelativistic 
degeneracy regimes due to increase in the plasma density. The variation in the potential-minimum location also follows 
quite dissimilar rout in the two degeneracy regimes. 

Figure 3 closely inspects the variations of Shukla-Eliasson attractive potential parameters, namely, the binding 
energy and binding length for the relativistic degeneracy regime relevant to the white dwarf core density for three 
different core compositions of Helium, Carbon and Magnesium. It is observed that the unlike for the binding energy, 
4>{fmin)i the binding length, f m i n , is unaffected by the change in plasma atomic-number, Z. It is further remarked 
that in the relativistic densities the binding length decreases suddenly with increase of the plasma density and reaches 
an asymptotic value of about I0 _14 cm, which is which is two orders of magnitude smaller that the nucleon diameter 
which may seem controversial. However, One should note that the hydrodynamic theory used in this model may not 
even be valid beyond the quantum screening length. On the other hand, other processes like nuclear fusion takes place 
before such short ionic distances are met. Figure 3(b) within a fixed relativistic density heavier atoms crystallize into 
stronger lattices. Crystallization in white dwarf cores are sometimes wrongly attributed to the presence of strong 
gravitational inward pressure at the center of the star. However, such confinement of ions in the core may only lead 
to quasi-lattice formation with viscous ion fluid due to lowering of ion coupling parameter value in fusion-medium 
temperatures. Slattery et al. 54j have shown that for ion coupling parameter values of T < 178 the plasma will be in 



liquid state. Therefore, Shukla-Eliasson attractive force can attribute appropriately to the perfect crystallization in a 
wide density range relevant to both cores of small planets as well as that of a typical white dwarf star. 

Figure 4 reveals yet another fundamental feature of Shukla-Eliasson screening for relativistic density regimes. 
It depicts variation of the normalized binding strength versus the normalized binding length for different plasma 
compositions. It is observed that with decrease of the binding length the binding energy increases for all plasma 
composition with more pronounced increase for higher plasma number density. Such variations reveal that the plasma 
crystals with relativistic density are unstable under external compression. Any force exerted from outside such as the 
gravitational inward pull for the case of a compact star would decrease the binding length and lead to collapse of the 
configuration. Such progressive increase in the plasma binding energy due to the gravitational pressure would also 
result in the energy release in the core leading the nuclear fusion ignition. As the atoms fuse they will form heavier 
elements which are lead to form even more stronger bounds and the release of more energy. Such loop-hole may 
proceed infinitely until the ultimate collapse of the white dwarf core releasing enough energy to lead to a cataclysmic 
supernova to take place. It should be noted that such mechanism called the spontaneous core collapse is only possible 
in the relativistic degeneracy regime, r s > r S 2(Z) and is not possible for r s < r s \{Z). 
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FIG. 2. (Color online) Variations of normalized electric potential (p — rsfylZe- for a > 1/4 with changes in plasma density 
and the ion charge-state for nonrelativistic (Figs, a and c) and relativistic (Figs, b and d) degeneracy density regimes. The 
distance f — r /rs is measured in the Bohr-radius unit. The thickness of curves appropriately reflect the variation in the varied 
plasma parameter, in each plot. 



We conclude that Shukla-Eliasson quantum screening potential which goes beyond the ordinary Debye or the 
Thomas-Fermi models can successfully explain the structure formation in solid-density matter, cores of small planets 
and compact stellar objects. Such theory which relies on the complete quantum interpretation of collective fermion 
behavior can even lead to deep understanding of the stellar formation and structure. It can also result in breakthroughs 
in pure as well rapidly growing nanotechnology sciences. 
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FIG. 3. (Color online) Variation in the binding energy and length in the relativistic degeneracy regime for three different 
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FIG. 4. (Color online) Variation of the normalized plasma binding energy in terms of the binding length for various plasma 
composition, namely, Helium, Carbon and Magnesium. The plot shows a feature specific to relativistic degeneracy regime 
called the spontaneous core collapse relevant to white dwarf stars. The thickness of curves being a measure of the value of the 
varied parameter. 
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